Mean-field approach to ferromagnetism in (III,Mn)V diluted magnetic 
semiconductors at low carrier densities 
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We present a detailed study, within the mean-field approximation, of an impurity band model for 
III-V diluted magnetic semiconductors. Such a model should be relevant at low carrier densities, 
below and near the metal-insulator transition. Positional disorder of the magnetic impurities inside 
the host semiconductor is shown to have observable consequences for the shape of the magnetization 
curve. Below the critical temperature the magnetization is spatially inhomogeneous, leading to very 
unusual temperature dependence of the average magnetization as well as specific heat. Disorder is 
also found to enhance the ferromagnetic transition temperature. Unusual spin and charge transport 
is implied. 
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I. INTRODUCTION 

Diluted magnetic semiconductors (DMS) are semicon- 
ductors of the general type Ai_a;Mj;B, where AB is ei- 
ther a II- VI or a III-V semiconductor and M a magnetic 
element, most commonly Mn. Substitution of a small 
fraction x of the element A by Mn impurities (and in the 
case of II- VI semiconductors an additional charge dopant, 
such as P on the B site) leads to the emergence of a semi- 
conductor with ferromagnetic properties.u This is due to 
the interactions of the S — ^ Mn spins (coming from the 
half-filled 3d shell of the Mn) with the spins of the charge 
carriers introduced by the Mn dopants, or, in the case of 
II- VI semiconductors, by the additional dopant. This 
opens up the possibility of manipulating (through dop- 
ing) not only the charge, but also the spin properties of 
the semiconductor. The ability to control the properties 
of a system by acting on the spin of its charge carriers 
is the subject of the new field of "spintronics" , and is 
believed to hold the promise to developing devices which 
combine storage functionalities (such as memory devices) 
together with information processing functionalities. 

The recent demonstration of Curie temperatures of the 
order of 100 K in Gai_a;Mna;As samples with x « 0.05, 
grown at low-temperatures using molecular-beam epi- 
taxy (MBE) techniques ,□□ has heightened the interest 
in understanding the physics of these alloys. By now, 
it is well established that the main magnetic interac- 
tion is an antiferromagnetic (AFM) exchange between 
the Mn spins and the charge carrier spins. As a re- 
sult, an effective ferromagnetic (FM) interaction arises 
between the Mn spins through carrier-induced ferromag- 
netism. Several theories, most notably Ruderman-Kittel- 
Kasuya-Yosida (RKKY),eI have been used to explain, 
within a mean-field approximation, this phenomenon. 
More recently, dynamic correlations as well as arbitrarji 
itinerant-carrier spin polarizations have been included.^ 
However, all these models assume that the charge carriers 
(the holes) occupy a Fermi sea in the valence band. This 
implicitly assumes a spatially homogeneous distribution 



of the holes throughout the entire system, and therefore 
completely neglects the role of the Coulomb attraction 
between the charge carriers and the Mn ions, as well as 
the role of disorder which is always present in such alloy 
systems. j-. 

It is well knowru that Mn impurities in a III-V semi- 
conductor create a trapping potential for holes. The as- 
sociated acceptor levels are about 100 meV above the 
top of the valence band, for GaAs, and interjactions will 
lead to their splitting into an impurity band.LHj At T = 0, 
the holes occupy the states of lowest energy and therefore 
they first occupy states in the impurity band. Only if the 
Fermi energy (or thermal energy ksT) is large enough are 
states in the valence band occupied as well. It has been 
found experimentally that these alloys are heavily com- 
pensated, leading to rather small hole concentrations, of 
the order of 10% of the Mn concentration.EJ The cor- 
responding small Fermi energy implies a long screening 
length for the Coulomb interactions, and opens up the 
possibility that holes are actually moving (through hop- 
ping processes) in the impurity band formed of states 
localized about the Mn impurities. Strong experimental 
evidence for this scenario comes from electrical conduc- 
tivity measurements, which at low tenrneratures reveal 
Mott variable-range hopping behavior. Q^tll As a result, 
randomness in the position of the Mn spins (and asso- 
ciated random potential) could be expected to play an 
important role. 

In this paper we investigate, at the mean-field level, a 
model in which holes move in a band formed of impurity 
states, neglecting the existence of the valence band states. 
Some results of this model have already been reported in 
Ref. 12. While a complete model should include both 
impurity and valence band states, we believe that our 
model is a good first approximation, especially at low 
temperatures, since we find typical Fermi energies of a 
few tens of meV, smaller than the average ~ lOOmeV 
between the impurity band and the valence band. 

The paper is organized as follows. In Section II we de- 
scribe the model Hamiltonian, the self-consistent mean- 
field approximation as well as the specific parameters 
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used in the numerical calculations. It is important to 
emphasize that although we chose numbers from the lit- 
erature specific to the GaMnAs problem for illustration 
purposes, similar arguments and physics would hold for 
other III-V compounds. In Section III we present our re- 
sults. We first analyze the hypothetical case in which all 
Mn impurities are ordered in a simple cubic lattice. This 
allows us to understand the unusual shape of the mag- 
netization curves we obtain. More importantly, it allows 
us to clearly identify the non-trivial effects of random- 
ness (disorder) in the Mn positions. These are analyzed 
in some detail and a clear physical picture of their im- 
portance to the problem emerges. We also analyze the 
metal-insulator transition, the effects of on-site random 
disorder, as well as other possible interactions. The sen- 
sitivity of our results to variations in the impurity band 
parameters is also discussed. Finally, Section IV contains 
a summary and conclusions. 

II. MODEL 

Based on experimental work in the liter ature,lil we as- 
sume that Mn only substitutes for the group III ele- 
ment of the III-V semiconductor. The III-V semicon- 
ductor is assumed to have a zinc-blende structure. Let 
Ri,i = 1, A^d be the random positions of the Nj, Mn 
dopants in the face centered cubic (FCC) sublattice of 
the group III element. Each Mn impurity is associated 
with a spin- 1 from its half-filled 3d shell. Since Mn has 
nominal valence II, when it substitutes for the valence-III 
element it will act as an acceptor. As a result, it can trap 
a hole shallow level, characterized by a hydrogenic Bohr 
radius as- Let Nh be the total number of holes trapped 
about various Mn sites, and p = Nh/Nd the relative hole 
concentration. A hole can hop from one Mn impurity to 
another one, while its spin is antifcrromagnctically cou- 
pled to the Mn spins in its vicinity. 

For simplicity, although the charge carriers of this sys- 
tem are holes, in the following we use an "electron" - 
formalism to analyze them. In other words, we in fact 
analyze an equivalent system doped with hypothetical 
donors, with impurity levels below a conduction-like 
band, instead of being above a valence-like band. The 
main difference between the two models is (i) the hole 
spin is I and (ii) the envelope wave function of hole has 
cubic symmetry, rather than spherical symmetry of a 
donor. This leads to some quantitative differences (as 
has been shown by MacDonald and coworkers for free 
holes in Ref. 13) but the essential aspects of the prob- 
lem, namely (a) disorder and (b) the effect of impurity 
potentials, which we concentrate on here, remain sub- 
stantially unaffected. 

The Hamiltonian we study is: 

i,3 i 
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Here, is the creation operator of a charge carrier with 
spin a in the bound orbital associated with the ith Mn 
impurity. The first term in Eq. (|^) describes charge 
carriers hopping between various Mn sites, with the 
hopping matrix tij dependent on the hopping distance 
r = \Ri — Rj\. We take an exponential form like for hjz^ 
drogcnic orbitals t{r) = 2 (1 + r/aB) exp {—r/aB) Ry.Efl 
We present detailed results for this case; however we dis- 
cuss in subsection F the changes that result with other 
models of the hopping integral. From these different cases 
we discern the universal aspects of the impurity models 
and also the requisite conditions for the existence of the 
experimentally observed ferromagnetism. 

The Rydberg (Ry) is defined by the binding energy of 
the charge carrier to the shallow trap, E},. The heav- 
ily compensated nature of the system gives rise to ran- 
dom potentials coming from the charged centers respon- 
sible for compensation. We model this through an on- 
site random potential u{i) which leads to the second 
term in Eq. |l|. A full model would require a self- 
consistent determination of this potential taking in con- 
sideration screening processes. The third term is the 
AFM interaction between Mn spins S{i) and charge car- 
rier spins. Since the Mn spins are very localized (the 
3d shell has a radius of the order 1-2A) the antiferro- 
magnetic exchange integral is taken to be of the form 

Jij — J|'/'(j,j)P = Jexp 2|^i — _Rj|/aB^ , refiecting 

the probability of finding the charge carrier in the s-type 
shallow level about impurity j near the ith Mn spin. In 
this notation, J is simply the AFM exchange integral of 
an isolated localized hole centered at a Mn impurity with 
the spin of the Mn. The last line in Eq. (^ describes the 
interaction of the carrier and Mn spins with an external 
magnetic field H. 

A. Mean-Field Approximation 

We treat the antiferromagnetic interaction at the 
mean-field level, given by the factorization 

S{i) ■ (J J ^ 5M„(*)fT| + S'{i)sh{]) - SMn{i)sh{j), (2) 

where 5j\/„(i) = {S^{i)) and Sh{j) = {crp are the ex- 
pectation values of the ith Mn spin and of the total spin 
created by charge carriers at the jth Mn site, respectively, 
and must be computed self-consistently at each site. For 
simplicity of notation, we used (jj = c^ja\^a(3Cji3 in Eq. 
(^). In writing Eq. (H) we made the implicit assumption 
that the rotational symmetry is broken in the direction 
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of the external magnetic field H (defining the z-axis) for 
all spins. In fact, one may start with a Heisenberg-like 
factorization in the absence of external magnetic fields, 
but we find that in the self-consistent configurations the 
spins are always coUinear, and as a result we regain the 
Ising-like factorization of Eq. (||) . Disorder induced non- 
collinear ground-states have been recently suggested in a 
model which assimaes that holes occupy a Fermi sea in 
the valence band.tJ 

The mean-field Hamiltonian can be separated into 
three parts 

'HmF = 'Hspin + 'Hcc + Ti-const- (3) 

The spin Hamiltonian may be rewritten as 



where 



Ti-spin — — ^ HiS^{i), 



(4) 



(5) 



is the effective magnetic field for the Mn spin at site 
i, including the charge carrier contributions. Then, the 
average Mn spin at site i is 



(6) 



where Bsix) = {S + ^) coth [(5 + ^)x] - \ coth § is the 
Brillouin function corresponding to the spin S = 5/2, 
and /3 = l/kBT. The contribution to the internal energy 
from the spin Hamiltonian equals 

f/spm(T)--^i/,5A/„(i). (7) 

i 

The charge carrier Hamiltonian can be rewritten as 



(8) 



where 



= ? "^11 l0(^i)l''5A/„(j) - giJiBH -h u{i) (9) 



is the effective on-site energy created by the Mn spins, 
the on-site disorder and the external magnetic field. A 
chemical potential /i has been added since we treat the 
charge carriers in the grand-canonical ensemble. 

The charge carrier Hamiltonian can be diagonalized to 
obtain: 



using the linear combinations 



(10) 
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where ipntrii) is the probability amplitude to find a charge 
carrier occupying level n with spin a in the shallow state 
centered at site i. The chemical potential is given by the 
condition 



(12) 



where f{Ena) = [exp (/3(i?„cr — A^)) + 1] is the Fermi 
distribution describing the occupation probability of the 
level (ncr). 

The expectation value of various charge carrier opera- 
tors can be computed in a straightforward way. In par- 
ticular, the average spin created by charge carriers at site 
i is 

= J E [\^n^i^)\^fiEn^-) - li'nm^fiEnl)] , (13) 



while the contribution of the charge carriers to the total 
internal energy is 



U,,iT)=J2EnafiEna)- 



(14) 



Finally, the third term in Eq. (^) contains the constant 
terms from the mean-field factorization 

T-Cconst ^ \4'ii,j)\'^SMn(i)sh{j) = Uconst{T). (15) 

As a result, the total internal energy of the system is 
given by 

U{T) = -gi^BHY^SMnit) + J2E,-,,f{E^,), (16) 



where, in the absence of external magnetic fields, all the 
contribution comes from Ucc{T). However, the charge 
carrier Hamiltonian Ticc contains the interaction with the 
average Mn spins incorporated in the eigenenergies En„ , 
so in fact their contribution is also included in this term. 
The specific heat of the entire system is given by 



Cv{T) 



dU{T) 
dT ■ 



Various other quantities of interest can be computed in 
a similar fashion. 

We solve the mean- field equations (H)-(|l^) using an it- 
erative algorithm. We start with a guess for the initial 
Sunii) configuration for each temperature T of interest 
(details about this are provided in section III.B). We nu- 
merically diagonalize the charge carrier Hamiltonian Eq. 
(^) and find the charge carrier eigenvalues, eigenfunc- 
tions as well as the chemical potential from Eq. (p^). 
This allows us to compute the expectation values for the 
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charge carrier spins at each site Sh{j) [from Eq. (13)] 
and therefore obtain the effective magnetic field at each 
Mn site Hi [Eq. ^] . The new expectation values for the 
Mn spins at each site SMn{i) are then obtained from Eq. 
(^ and the iterations are repeated until self-consistency 
is achieved. We define self-consistency as corresponding 
to the situation where the largest absolute variation of 
the on-site charge carrier energies eio- [Eq. (||)] between 
two successive iterations is less than 10"'^. This corre- 
sponds to relative errors of 10~^ or less for all computed 
quantities. 



B. Parameters 

We consider a NxNxN FCC sublattice of the valence- 
III element, of lattice constant a (a = 5. 65 A for GaAs). 
We assume throughout this paper that the system has 
periodic boundary conditions. The Mn doping is charac- 
terized hy X = Nd/4:N^, leading to a Mn concentration 
CMn = 'ix/a^. The hole concentration is Ch = pcMn, 
where p = A^^fi/TV^i-^/alues of interest are x — 0.01 — 0.05 
and p = 5 — 10%.li2l The choice of the system size N is 
described in the following subsection. 

Throughout this paper we use parameters specific to 
the Gai_2;Mna;As system. The binding energy of the hole 
(defining the Ry unit of this problem) is Eh = 112.4 
meV=lRy.l] Using the Jjuttinger Hamiltonian in the 
spherical approximationQ we find the effective mass for 
the heavy hole to be 



71 - (673 + 472)75 



0.56 nip 



(17) 



where the 7-coefficients for GaAs are 71 = 7.65, 72 = 2.41 
and 73 = 3.28.0 This allows us to estimate the Bohr 
radius of the isolated impurity state aa3 



as 



y/{2mhEb) 



7.8 A. 



(18) 



This is in excellent agreement with another possible esti- 
mate, obtained by assuming that the hole is bound to 
its impurity by a pure Coulomb attraction, in which 
case Eb = e^/(2eas). Using the value e — 10.66 for 
GaAs leads to = 7. 82 A. While this agreement is 
probably jfortuitous, similar values have been used in 
literaturefl The characteristic value of the hopping inte- 
gral is t{4aB) = 20meV. Besides the binding energy Ei,, a 
second energy scale is provided by the exchange integral. 
We use the value J = 3e = 15 meV, where e = 5 meV 
is the value obtained in Ref. 7 for the antiferromagnetic 
interaction of an isolated hole with the spin of its own 
trapping Mn impurity. We include the factor of 3 as the 
simplest way to account for the fact that the heavy holes 
have spin projections jz = ±|, while in our model they 
are modeled as = ±i objects. The final parameter. 



W, (or its dimensionless counterpart W/Ei,) character- 
izes the on-site disorder due to uncompensated impuri- 
ties. We assume that u{i) has a uniform distribution in 
a range [— V7, VF]. Following Ref. 8 we find an estimate 
of W given by e^/ef, where f ^ l/^-i/n roughly the 
average distance between Mn impurities. For the typical 
charge carrier concentration considered Ch ~ 1.5 x 10^*^ 
cm"'^ and p = 10%, we find f 9 A ^ ub, which sug- 
gests W ^ 1 Ry. As we show later, while the magnitude 
of on-site disorder W influences the shape of the mag- 
netization curve, the critical temperature has only a 
very weak dependence on it. As a result, and given the 
fact that the compensation mechanism in these systems 
has not yet been fully clarified, we do not attempt a more 
detailed modeling of the on-site disorder at the present 
stage. 

Thus, in the absence of external magnetic fields, the 
problem depends on five dimensionless parameters, J/i?&, 
as/a, n/ifl'g, x and W/Ei. 



III. RESULTS 

A. Simple cubic superlattice of Mn 

In order to gain some insight in the behavior of the 
system, we first consider ths simplified case of Mn impu- 
rities placed in a simple cubic structure, with a super- 
lattice constant = a/{AxY^^. Strictly speaking, only 
concentrations x for which is commensurate with a 
would be physically acceptable. However, since we per- 
form this calculation only to get a feeling for the homoge- 
neous solution, we disregard the underlying GaAs lattice 
in this particular case and assume that Mn spins could 
be located anywhere in space. We also set the on-site dis- 
order u{i) = 0, and turn off the external magnetic field 
(i/ = 0). 

For the ordered case, translational symmetry implies 
Sh{i) = Stole, SMn{i) = Smv, for all sites i of the cubic 
Mn superlattice. As a result, the charge carrier Hamilto- 
nian Tiec is diagonalized with plane waves, and the self- 
consistent equations (p|)-(p^) reduce to: 



Smt. 



Shole 
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-B{PJeffShole), 



where the chemical potential is determined from 
Nh 1 



P 



7 dkT.f(Ek<r) 



Nd 87r3 

and the charge carrier eigenenergies are given by 
Ej:^ = e{k) + ^JeffSnn- 



(19) 



(20) 



(21) 



(22) 
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Here, k is measured in units of I/ol, and the integrals in 
Eqs. (^ ) and (|2l|) are performed over the first Brillouin 
zone — TT < ka < TT, a = x,y, z. The non-interacting 
charge carrier dispersion relation is given by 

e(fc)=^t(r,)e-^'-^^ (23) 

and 

Je// = J^|(/.(j)|2 = J^e"^. (24) 

i i 

In Eq. (p3| ) and ( p^ ) the sums are performed over the 
whole crystal, but the exponential decay of | </<(«) P and 
t{ri) lead to finite results. 

The self-consistent values of < Smu < 2.5 and the 
average charge carrier spin —0.5 < Sh = Sh{i)/Nh — 
ShoieNd/Nh < obtained for x = 0.01-0.05 andp = 10% 
are shown in Fig.|^. The overall signs indicate the AFM 
alignment of the charge carrier and Mn spins. 

The total magnetization of the sample, obtained by 
adding the Mn and charge carrier contributions, looks 
similar to the Smu curve, since the number of Mn spins 
is 1/p ~ 10 times larger than the number of charge car- 
rier spins, and they also have larger g- factors. Thus, we 
see that the magnetization curve does not have the typ- 
ical form of ferromagnetic systems. In fact, each curve 
shows three different regimes. Below Tc there is a re- 
gion where neither the charge carrier nor the Mn spins 
are yet fully polarized. The gap JeffSMn between the 
a =1 and a =t charge carrier bands (see Eq. (p2|)) in- 
creases quickly as the temperature decreases. This leads 
to a polarization of the charge carriers, which in turn 
polarizes the Mn spins even more. Since the number of 
charge carriers is relatively small, they are the first to 
fully polarize, at a characteristic temperature defined by 
JeffSMn — Ep ~ ksT. Here, Ep is the Fermi energy of 
charge carriers measured from the bottom of the a =[ 
band, and the condition simply means that the gap from 
the highest occupied a =[ level to the first available a =t 
level is larger than the thermal energy. As a result, below 
this temperature charge carriers are fully spin-polarized, 
Shoie = psh = —O.bp. From Eq. ( pO| ) we see that be- 
low this temperature, the Mn spins behave as if they 
are in a constant external magnetic field of magnitude 
H = JeffShoie- For temperatures such that PH <C 1, 
the Brillouin function may be linearized and we find that 
Smu ~ Jeffs/ksT (the Curie law), and Smu increases 
roughly like 1/T as the temperature decreases. This ex- 
plains the uncharacteristic concave upward shape of the 
Mn spin magnetization. Finally, below temperatures for 

which BsiPJeffs) « S (i.e. kgT < SpY^^e'^ meV), 
the Mn spins are also fully polarized. 

In the inset of Figj^ we plot the specific heat per Mn 
impurity for the same parameters. In all cases we see two 
distinct contributions. The lower peak is entirely due to 
the Mn spins, while the upper one is the charge carrier 



contribution. At low temperatures the charge carriers 
are all "frozen" at the bottom of the a =| band, and all 
the entropy is due to fluctuations in the Mn spins. This 
can be easily checked by computing (7^™'* = dUspin/dT, 
with Sfioie = —0.5p substituted in Eq. (j7|). This accounts 
for the entire lower peak. At the higher temperatures 
the Mn spins are almost free (the effective magnetic field 
orienting them is very small), and therefore right below 
Tc the entropy is dominated by spin fluctuations of the 
charge carriers. 

For increasing charge carrier concentrations p, one ex- 
pects that the temperature where charge carriers become 
fully polarized decreases (since Ep increases) and thus 
the unusual regime with fully polarized charge carriers 
and Smu ^/ksT is restricted to smaller intervals. In 
other words, we expect that for larger p values the Mn 
spin magnetization should begin to look more like the 
characteristic convex upward (concave downward) form 
seen in usual ferromagnetic systems. This is confirmed 
in Fig. ^, where the average charge carrier and Mn spins 
are plotted for x = 0.02 and p = 5, 10, 25 and 40%. 

It is interesting to note that the critical temperatures 
obtained for this homogeneous case using the nominal 
parameters and Bohr radii from the literature are ac- 
tually in good agreement with experimentally measured 
values. In Fig. y we show the critical temperatures for 
three GaMnAs samplesjl2l which are estimated to have 
p = 5 — 10%. The experimental points fall right in be- 
tween the theoretical curves corresponding to the two 
charge carrier concentrations p. However, it is important 
to emphasize the fact that fairly small variations in any 
of the parameters can lead to rather large variations in 
Tc. Indeed, from Fig. ^ we see how sensitive Tc and 
the shape of the magnetization are to p. If we increase 
the Bohr radius by just lA, while keeping all the other 
parameters fixed, the critical temperatures increase by 
roughly 50%, and the experimental points are well below 
the new p — 5% theoretical estimations. In the following 
section we show that disorder in Mn positions has, at 
least at the mean- field level, a large effect on Tc- At the 
same time, the mean-field approximation itself underes- 
timates thermal fluctuations, and therefore may substan- 
tially overestimate critical temperatures. Finally, the ac- 
tual transition temperature varies substantially with the 
precise form for the hopping parameter (see section F). 
Therefore, we conclude that the good agreement shown 
in Fig. ^ is likely fortuitous. 

B. Effects of disorder in Mn positions 

We now analyze the effects of randomness in the Mn 
impurity positions on the shape of the magnetization 
curve and the value of the critical temperature. We 
restrict ourselves in this section to the case of no on- 
site disorder, i.e. u{i) = 0, and vanishing external mag- 
netic field. In this case, the system is no longer homoge- 
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neous. As a result, we have to limit ourselves to a finite 
N X N X N FCC lattice with periodic boundary condi- 
tions, choose randomly the positions of the Mn spins 
and solve Eqs.(|5|)-(p^ self-consistently for each site. 

A technical question is how large a system wc should 
consider in order to assume that the "bulk limit" has 
been reached. The size N of the system is chosen so as to 
minimize finite size effects. These are monitored through 
both their effect on the magnetization curves (especially 
on Tc), as well as on the total density of states (DOS). We 
find that for x = 0.05, finite-size effects become negligi- 
ble for systems which have more than Nh ~ 50 holes. As 
a result, the minimum size we use for this concentration 
is iV = 14, corresponding to Nd — 548, = 55. For the 
smallest Mn concentration we investigated, x — 0.00926, 
we find that finite size effects are small for systems having 
as few as 12 holes (corresponding to = 15, close to the 
previous value). However, even for this lower concentra- 
tion, most of the results shown are obtained for lattices 
of size N — 24, corresponding to Nd — 512, — 51. 
This ensures that grand-canonical fluctuations in the to- 
tal number of charge carriers are minimized as well. 

To solve the equations ([^)-([l3|) self-consistently, we 
start with different initial Mn spin values SMn{i) to use in 
the charge carrier Hamiltonian Ticc for the first iteration 
of the process to self-consistency. We first consider biased 
initial conditions. In this case, we start the first iteration 
for the lowest temperature considered by assuming that 
all Mn spins are fully polarized, SMn{i) = 2.5. After 
several iterations, the self-consistent values SMn{i) cor- 
responding to this temperature are found. We then use 
these values as the initial guess for the next higher tem- 
perature considered, etc. This allows us to find, at each 
temperature, the self-consistent solution with the highest 
possible total magnetization. We next start each search 
with random values for both the magnitude and the sign 
of SMn{i) for each temperature considered. Since in 
principle there could be metastable solutions, the "true" 
mean-field solution is the one corresponding to the lowest 
total free energy. 

In Fig.^, left panel, we show the expectation values 
for the average Mn spin Smu = X^i Sunii) and the 
average charge carrier spin Sft, = 1 /Nh (i) obtained 

for one disorder realization using biased (full lines) and 
random initial configurations (circles) for p = 10% and 
X = 0.00926. For random initial conditions we actually 
plot \SMn\ > (full circles) and — |s;i| < (empty cir- 
cles), since the two expectation values always have op- 
posite sign, but the orientation is arbitrary for random 
initial conditions. Random initial conditions lead to var- 
ious self-consistent configurations, with magnetizations 
smaller or equal to the maximum possible value given by 
the biased configuration. The smaller value of average 
magnetizations for the random initial conditions config- 
urations is not a consequence of smaller polarizations of 
individual spins, but of the appearance of regions with 
local magnetizations pointing in different directions. The 
right panel of Fig.0 shows the difference between the total 



energy per Mn spin obtained with random initial condi- 
tions, and that of the biased configurations. At low tem- 
peratures, all the random configurations are much higher 
in energy than the biased configuration, with the differ- 
ence increasing for configurations with lower total magne- 
tization. However, by ksT/J ^ 0.5 the energy difference 
between configurations becomes comparable to ksT, as 
evident in the intersection with the solid line which is a 
plot of ksT/Nd- At these temperatures thermal fluctua- 
tions will enable the system to vary continuously among 
these various states, effectively suppressing the magneti- 
zation and therefore lowering T,.. A proper treatment of 
the effect of thermal fluctuations requires going beyond 
the mean-field approximation, e.g. by a Monte Carlo 
simulation. 

We have found similar results for various other Mn and 
hole densities. We have looked in all at over 150 samples 
of varying sizes N — 14 — 24 corresponding to Mn con- 
centrations X = 0.00926 — 0.05 and relative hole concen- 
trations p ~ 10 — 30%. In all cases, the most polarized 
(biased) state has been found to have the lowest total free 
energy in our model. For larger x, we in fact find that 
most random initial samples converge to the biased limit, 
signifying a more robust magnetization than at lower x. 
We conclude that for this range of concentrations the 
system is indeed ferromagnetic at low temperatures, and 
that the biased configuration curves may be used to ob- 
tain the lowest energy mean-field configuration possible. 
However, the Tc given by this mean-field biased curve is 
likely to be significantly larger than the one provided by 
a Monte-Carlo simulation, with the difference likely to 
grow as the concentration decreases. 

Comparing typical Tc values obtained for random Mn 
configurations (Fig. ^) with those obtained for the or- 
dered Mn lattices with similar x and p values (see Fig. 
|l|) we see that the shape of the magnetization curves is 
significantly changed by randomness. Tc is increased, 
while the curves become even more concave. In fact, the 
Mn spins do not reach the saturation limit S'a/„ = 2.5 
until very low temperatures. While at first sight this sig- 
nificant increase of Tc for the disordered case may seem 
puzzling, in fact it has a very physical explanation. In 
the disordered sample there are regions of higher local 
concentration of Mn. The charge carriers prefer these re- 
gions, since they can lower their kinetic energy by mov- 
ing among several nearby Mn sites. They can also lower 
their magnetic energy by polarizing the spins of these 
Mn impurities. As a result, these regions of higher Mn 
concentration will become spin-polarized at higher tem- 
peratures than the average sample, pushing Tc up. 

This point can be illustrated by looking at a histogram 
of the total density of charge carriers at site i, defined as 

pit) = \Hhj)\'{c]^cn + c.Vji)- (25) 
j 

[In other words, p(i) is given by the probability p{i) = 
{cltCii+clfCii) of finding a charge carrier at site i, plus ex- 
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ponentially small contributions due to tailing from charge 
carriers found on nearby Mn sites]. Since at low tem- 
peratures only spin down charge carrier states are sig- 
nificantly occupied, we can also interpret p(i) as being 
proportional to the effective magnetic field acting on the 
ith Mn spins [see Eq. (||)], with external magnetic field 
H = 0). In Fig. |§i we show, on a logarithmic plot, such 
a histogram obtained for 25 random Mn configurations 
with X = 0.00926 and p = 10% at ksT/J = 0.0006 (dot- 
ted line) and ksT/J = 0.6 (full line). The vertical line 
indicates the position of the i5-function for the histogram 
corresponding to the homogeneous system (ordered Mn 
superlattice). In that case the density p{i) at every site 
is the same, and is given by p = p exp (— 2ri/aB). 
For small concentrations x the sum is roughly equal to 
unity, and the charge carrier density and effective mag- 
netic field at each Mn site is p, respectively Jp. On the 
other hand, for the random (disordered) configurations, 
a double-peaked structure is clearly visible. There is one 
sharp peak centered about p ~ 0.6 > p corresponding 
to densities much higher than the average, and a second 
much broader peak centered at exponentially small values 
p « 10~^ <C p. In other words, there are some Mn sites 
with a large charge carrier density, which strongly po- 
larizes the respective Mn spins up to high temperatures. 
These correspond to the high density regions of Mn, and 
define a much higher than the one of the homogeneous 
system. The rest of the Mn sites have a small charge 
carrier density (or effective magnetic field) and therefore 
they very quickly become depolarized as the temperature 
increases, leading to the fast decrease in the average Mn 
spin value Smu- This phenomenology can be captured 
quite accurately by dividing the spins into strongly and 
weakly interacting ones, depending on whether their ef- 
fective magnetic field is larger or smaller than the corre- 
sponding thermal energy.E3 A mix of ferromagnetic and 
paramagnetic contributions to the M{H) curves, which 
can be attributed to strongly, respectively weakly mter- 
acting spins, has also been observed experimentally.EII As 
the temperature increases to fc^T/J = 0.6 (just below 
for this density), the double-peaked structure is still ap- 
parent, although it becomes more centered around the 
average value and the peak corresponding to strongly in- 
teracting spins decreases. For temperatures well above 
Tc (see inset of Fig. ||a) the distribution width decreases 
even more, although it still extends over more than two 
orders of magnitude. This behavior suggests that as the 
temperature is lowered through Tc, the charge carriers 
start to polarize the most dense clusters of Mn. As a re- 
sult, their wave- functions become more concentrated in 
these high-density Mn areas, where the carriers can fur- 
ther lower their magnetic exchange energy. This leads 
to the increased weight of the higher p{i) peak. The 
concentration of charge carriers in the high-density ar- 
eas implies a further depopulation of the low Mn density 
areas, pushing the low edge of the p{i) distribution to- 
wards lower values. With decreasing temperature the his- 
togram quickly changes and for ksT/J < 0.3 it already 



has a shape identical to the that of the low temperature 
ksT/J = 0.0006 histogram. 

In Fig. Ijb we compare the low temperature histograms 
for different concentrations x = 0.00926 (dotted lines) 
and X — 0.05 (full lines). The vertical lines again show 
the corresponding values for the homogeneous (ordered) 
systems. For x = 0.05 the double-peak structure is also 
clearly visible. However, the whole histogram shifts to 
higher p{i) values. This is consistent with the fact that 
for larger charge carrier concentrations the overall inter- 
actions are increased. [It is interesting that for x — 0.05 
there is a finite concentration of sites for which p{i) > 1 
(logio[p(«)] > 0). Since the probability of finding a hole 
at any site p{i) = (c-|Ci| + c-|Cix) < 1, this suggests that 
in this case some Mn spins strongly interact with sev- 
eral charge carriers that are nearby [see Eq. (p5|)]. For 
X = 0.00926, however, p{i) < 1 for all sites, suggesting 
that Mn spins interact at most with one charge carrier]. 

To further check this picture of enhancement of Tc 
within our model, we have "tuned" the amount of disor- 
der (randomness) in the Mn positions. Fig. ^ shows 
curves of average charge carrier and Mn spins as a 
function of temperature for four different distributions 
of Mn spins. They all correspond to the same values 
X = 0.00926 and p = 10%. The curve with the low- 
est Tc is the curve for the ordered Mn superlattice, with 
gl = 3a. The next curve corresponds to a weak disor- 
der configuration in which each Mn atom is allowed to 
randomly choose one of the 12 nearest neighbors of the 
underlying FCC sublatticc. In other words some random- 
ness has been allowed for, although the configuration is 
still quite homogeneous, with one Mn site inside each 
cubic supercell. Even this small amount of disorder is 
seen to have a significant effect on the shape of the mag- 
netization and the mean-field Tj. value. Some Mn spins 
now have a higher effective charge carrier concentration, 
pushing Tc to higher values. However, at low tempera- 
tures we see that the average Mn spin is smaller than that 
of the ordered lattice, meaning that some Mn spins are 
in much lower effective magnetic fields (have lower local 
charge carrier density) and only get saturated at much 
lower temperatures. The third curve corresponds to a 
medium disorder configuration in which Mn impurities 
are allowed to pick any sites on the FCC sublattice, as 
long as the distance between any two of them is at least 
2a. This allows for even more randomness, leading to 
even higher Tc, while at low temperatures the average Mn 
spin is even more suppressed. Finally, the curve with the 
highest Tc corresponds to a completely random (strong 
disorder) Mn configuration. Monitoring the histogram 
of on-site densities for these cases, we find the expected 
behavior: increasing disorder leads to a larger spread of 
the densities about the average value of the ordered su- 
perlattice, leading to both the increase of mean-field Tc 
as well as the decrease of the saturation temperature. 

A qualitatively similar picture holds for higher Mn con- 
centrations as well as higher hole densities, although the 
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effects are quantitatively substantially less. In Fig. |^ we 
show the Mn and hole spins for both the simple cubic su- 
perlattice and the random Mn distribution on the FCC 
sublattice for x = 0.05 for two different p. While Tc 
is again substantially larger in the random system, the 
percentage increase is smaller than in the x — 0.00926 
case. Increasing the hole concentration from p — 10% 
to p = 30% makes the curve shape more conventional 
(convex upward). The reason is simply that the fluctua- 
tions in the local doping are smaller at higher Mn concen- 
trations, and increased hole doping further reduces the 
width of the density distribution. 

All the curves shown so far correspond to one partic- 
ular random distribution of the Mn impurities on the 
FCC sublattice. As the amount of disorder in the po- 
sitions of the Mn sites increases, so do the variations 
between curves corresponding to various realizations of 
disorder in this mean-field approach. In Fig. H we show 
a typical spread for 25 different disorder realizations for 
both X = 0.00926 (upper panel) and x — 0.05 (lower 
panel). In both cases we see that while the curves have 
similar shapes, there is a significant variation present. In 
particular, for the lower concentration, we see that most 
magnetization curves have a very elongated tail near Tc, 
arising from rare dense clusters of (usually nearest neigh- 
bor) Mn spins that are polarized by the holes. This tail 
will likely be destroyed by thermal fluctuations if one 
goes beyond the mean-field treatment. For the higher 
concentration most curves have ksTc ^ 0.9J, although 
again there are a few which have longer tails arising from 
magnetization of dense clusters of Mn impurities. 



C. Metal-Insulator Transition 

According to Mott's criterion, a doped semiconduc- 
tor goes through a metal-insulator transition (MIT) for 

a charge carrier concentration given by n\l^aB ~ 0.25. 
(For compensated systems, the critical density is ex- 
perimentally found to be somewhat larger). Neglect- 
ing the effects of compensation and assuming p = 10%, 
for as = 7. 8 A this corresponds to a Mn concentration 
x ~ 0.015. Given the variations in the compensation pa- 
rameter p for different samples, this is in good agreement 
with experiimental measurements which indicate a MIT 
at x - 0.03.N 

One way of monitoring the MIT is by determining 
whether the charge carrier states in the vicinity of the 
Fermi level are localized or extended. We characterize 
the charge carrier states using the Inverse Participation 
Ratio (IPR), defined for each state (ncr) by 



IPR{na) 



For a state extended over the Nd sites of the system, one 
expects |V'ncr(i)| ~ l/\/Nd, and therefore IPR{na) ^ 



l/Nd- In other words, for extended states IPR{na) is 
inversely proportional to the size of the system. For lo- 
calized states, IPR{na) is inversely proportional to the 
number of sites over which the wavefunction is localized, 
and therefore independent of the size of the system. 

In Fig. H we plot IPR{E, a) for a completely disor- 
dered system with x = 0.00926 and p — 10%, at a low 
temperature {ksT/ J — 0.0006). For each system size we 
consider 100 different disorder realizations, and we aver- 
age the IPR of all the states with eigenvalues within 2 
meV of each other. We show both the cr =[ and <t =t 
sub-bands for three system sizes, Nd = 125, 512 and 1000, 
corresponding to Nh = 12, 51 and 100. A large gap of size 
4.2J = 63 meV opens between the two subbands, and 
only states in the a —[ are occupied at this low temper- 
ature. The position of the Fermi energy is shown by the 
vertical line, at about 45 meV above the bottom of the 
band. We clearly see that states at the bottom of either 
band are localized, with the IPR independent of the sys- 
tem size. At higher energy, however, the states become 
extended, with the IPR ~ l/Nd- The mobility edge for 
the (T —i band is at an energy of about -100 meV. The 
Fermi level is below the mobility edge, signifying that 
this system is still insulating, in agreement with exper- 
imental measurements. However, the IPR of the states 
at the Fermi energy is less than a factor of 2 larger than 
the IPR at the mobility edge for our sizes. This suggests 
that significant tunneling may occur in between various 
regions occupied by the holes, leading to alignment of po- 
larization of all the high-density regions. The inset shows 
the density of states for the two spin-polarized subbands. 
The curves corresponding to the three system sizes fall on 
top of each other, proving that the "bulk limit" is already 
reached. In all cases we investigated, the DOS has an ex- 
tremely long upper tail, only part of which is shown. Due 
to the small relative concentration of holes p — 10%, only 
states very close to the bottom of the impurity band are 
occupied. As already mentioned, this supports our as- 
sumption that neglecting band states (which are roughly 
110 meV above the impurity band) is a good starting 
approximation. 

In Fig. |l^ we show the IPR for the a =| subbands of 
configurations with medium disorder (satisfying the re- 
striction that the distance between any two impurity sites 
is larger than 2a). Again, 100 configurations each for sys- 
tem of size Nd = 125, 512 and 1000 are averaged. The 
Fermi energy of the medium-disorder system is shown by 
the dashed line. The (J =1 sub-band (which is completely 
empty) is not shown. For the medium-disorder system 
the IPR values are lower than those corresponding to the 
strongly disordered system. This is the expected behav- 
ior, since in the limit of no disorder all wave functions 
must become fully extended and the system is metallic. 
In fact, we observe that even for the medium disorder 
case the Fermi energy is just above the mobility edge. 

The histograms presented in Fig. |^a show that as the 
temperature increases the distribution of charge carriers 
in the system becomes somewhat more homogeneous (the 
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width of the distributions decreases). This suggests that 
the charge-carrier wave- functions become more extended 
at higher temperatures, and therefore the system is more 
"metalhc" . This is p. quahtative agreement with resis- 
tivity measuremcntala which show a larger resistivity at 
low temperatures than above for all low-density sam- 
ples. We conclude that while the exact value of p and 
nature of disorder in Mn positions play a crucial role, 
the X = 0.00926 sample is most likely insulating. 

IPR curves for x = 0.05 and completely random sam- 
ples (strong disorder) are shown in Fig. |ll|. In this case, 
we have used larger size systems with Nd = 548, 981 and 
1600 in order to avoid finite size effects. Again, we see 
that states at the bottom of either subband are local- 
ized, while states at higher energies are extended (their 
IPR scales with 1/Nd)- In this case, the system is clearly 
above the metal-insulator transition, in qualitative agree- 
ment with experiment. _. .. . 

Experimentally it has also been observedoEJIlj that 
samples with even higher concentrations {x > 0.07) be- 
come insulating again. However, these samples also seem 
to have a much lower relative charge carrier density p 
(see Ref. 18). Since the x = 0.05 system is just above 
the MIT, it is reasonable to assume that a significant 
decrease in the density of charge carriers (leading to a 
significant decrease in the Fermi energy) could move the 
Fermi level below the mobility edge and therefore be re- 
sponsible for the re-entrant insulating state. 

The fact that these systems are either metallic, or 
not too far from the MIT, is important for obtaining 
a large critical temperature. Delocalization of the elec- 
tronic wavefunction over several sites, which allows the 
Mn spins to effectively communicate with each other, is 
the essential ingredient which leads to alignment of the 
polarization of all the high density regions. The charge 
carriers hopping between various high-density areas will 
force the alignment of Mn spins in each region to be the 
same, in order to lower their kinetic energy. However, 
maximizing the critical temperature seems to require a 
fine balance: increasing the disorder leads to increased 
Tc, but also to increased localization. If the charge car- 
rier states become so localized that there is no tunneling 
in between high density occupied regions, than the direc- 
tion of polarization of each such region is uncorrelated 
with the direction of polarization of the other regions, 
and the average magnetization of the sample will vanish 
(Tc — !• 0). On the other hand, a very homogeneous sam- 
ple has extended charge carrier states, but Tc is lower 
since in such case all the Mn spins are in a similar "av- 
erage" environment. 

D. Effects of on-side disorder 

We now consider the role of the on-site disorder u(i) ^ 
in our model. While the nature of the heavy com- 
pensation is not elucidated at this point, one may as- 



sume that compensation is due to the annihilation of the 
holes by some type of defect, leading to appearance of 
a background of charged impurities. These would create 
an electric potential u{i) at all Mn sites. The simplest 
way to describe it is to assume that the on-site energies 
u{i) are distributed with equal probability in the interval 
[—W, W] . In Fig. |l^ we compare the average Mn and spin 
charge carrier magnetizations obtained from a random 
Mn configuration with x = 0.05 and p — 10%, for vari- 
ous values of the on-site disorder cut-off W/Et — 0, 0.5 
and 1. On-site disorder does not affect Tc considerably 
(in all the simulations we performed, we find that Tc 
decreases slowly with increasing W). However, on-site 
disorder changes the shape of the magnetization curve. 
Due to on-site disorder, some of the charge carriers from 
the high density regions are pushed away into the less 
populated regions. As a result, the magnetization near 
Tc (which is dominated by contributions from Mn in the 
high-density regions) is suppressed. On the other hand, 
the low-temperature magnetization, which is dominated 
by the spins in the low density regions, is increased ac- 
cordingly. It is interesting to notice that the magnetiza- 
tion of the sample with W=l Ry varies almost linearly 
with temperature. Such unustial, M(T) dependence has 
been observed experimentally.Qll3 

We have also investigated a more detailed model for 
on-site disorder, which assumes that compensation is 
entirely due to As antisites. When a valence-V As 
atom substitute for a valence-III Ga atom, its two ex- 
tra electrons effectively remove two holes from the im- 
purity band. If all the compensation is due to such pro- 
cesses, then the number of As antisites must be given 
by Nas = — p)/2. Each such As impurity has an 

effective charge +2e, and therefore will contribute an on- 
site Coulomb potential -f2e^/er at a Mn impurity site 
which is at a distance r from it. However, since the Mn 
ions also have effective ionic charge — e, the As poten- 
tial is screened (partially compensated) by the potential 
of the Mn impurities nearby it. One could use a de- 
tailed formula u{i) = - 1]^ e^/erj,Mn + J2k'^^'^ /^'^k,As, 
with the first term describing the Mn contribution and 
the second one describing the As antisite contribution. 
An alternative, simpler form, is to assume that each As 
antisite only contributes to the on-site potential u{i) of 
its two nearest Mn neighbors, with the contribution to 
the other Mn sites being screened out by the contribu- 
tion of these two nearest Mn sites. The two formulations 
are qualitatively equivalent. Given the absence of more 
detailed information about the exact nature of compen- 
sation and screening, we investigated the simpler model. 
In this case, after we randomly choose the locations for 
the Mn impurities, we select random positions on the Ga 
sublattice for the As antisites as well. We find the two 
closest Mn neighbors for each As antisite (with each Mn 
selected as neighbor only for its closest As antisite), com- 
pute the corresponding values for u(i) and then proceed 
with the calculation as described previously. Typically, 
the on-site interaction in this model leads to a substan- 
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tial decrease in (see Fig. |^). Also, the shape of 
magnetization curves becomes even more concave, with 
the larger change for the hole magnetization, which no 
longer reaches full polarization in the T = limit. 

E. Effects of external magnetic fields 

Finally, we consider the effect of an external magnetic 
field, in the absence of on-site interactions u{i) = 0. We 
assume that g = 2 for the Mn spins. The precise value of 
the g-factor for the holes is not important, since we find 
that the magnetization is not changed if we vary g within 
a reasonable range. This is a consequence of the fact that 
each hole strongly interacts with many Mn spins, and the 
external magnetic field is just a small perturbation to the 
effective on-site energy Cia experienced by holes [see Eq. 
(^]. On the other hand, the external magnetic field leads 
to a significant change in the effective magnetic field Hi 
[see Eq. (||)] of each Mn spin, since any Mn interacts with 
very few holes (or practically none, for weakly interacting 
Mn). In Fig. HJ we plot the magnetization of a disor- 
dered configuration with x — 0.05 and p = 10%, in the 
presence of an external magnetic field H = 0,5 and lOT. 
The external magnetic field leads to a significant increase 
of the Mn spin magnetization at all temperatures, since 
it polarizes the many weakly interacting spins. It also 
leads to a saturation of the magnetization for tempera- 
tures ksT < gfisH, as expected. The magnetization of 
the charge carriers is also increased (in magnitude) in the 
presence of the magnetic field. This may seem puzzling 
at first, since one would expect that the external mag- 
netic field would favor a fiip of the charge carrier spin 
from cr =J, to (7 —1, leading to a decrease of the charge 
carrier magnetization. However, this is again due to the 
fact that each hole interacts with many Mn spins. As the 
magnetization SMnii) of each Mn spin is increased by the 
magnetic field, the effective negative magnetic field felt 
by the holes increases, more than compensating the pos- 
itive external magnetic field H [see Eq. (1)]. Addition of 
random on-site disorder changes the shape of the magne- 
tization curves, in a similar fashion to the one presented 
in Fig.0 

F. Beyond the hydrogenic model 

As emphasized in the introduction, the impurity band 
model with hydrogen-like, exponentially decaying hop- 
ping parameters, leaves out many aspects of the true 
Hamiltonian in a system like Gai_a;Mna;As, especially 
with large As antisite defects, or similar causes for the 
large observed compensation. First, and probably fore- 
most, is that the Mn dopant is not exactly "shallow", 
with a binding energy of over 100 meV. This suggests 
that the true wavefunction is substantially affected by 



central cell corrections and spin effects, which could af- 
fect the hopping integrals considerably. Secondly, the 
magnitude of t(r) is based on a two-center formalism for 
spherically symmetric wavefunctions valid at intermedi- 
ate to large separations. Besides the obvious complica- 
tion of anisotropy of the true hole wavefunction, these 
could lead to substantial renormalization of the effective 
t(r) at the densities of interest, especially at the upper 
end few percent). More microscopic calculationsQ sug- 
gest a significant renormalization of the energies within 
the impurity band, compared to the simple two-center 
tight-binding picture. Consequently, we discuss the ef- 
fects of changing the hopping parameter t(r) in some de- 
tail below. The effect of other approximations made in 
our study, namely, the neglect of the random potential 
due to the compensating centers, carrier-carrier interac- 
tions, and the valence band states, are discussed following 
that. 

The effect of reducing the magnitude of the hopping 
can be achieved by changing the Bohr radius, or the pref- 
actor. A change in Bohr radius is merely a renormaliza- 
tion of the effective density, while the effect of changing 
the prefactor is tantamount to changing the exchange 
coupling in the opposite direction, and then rescaling 
the temperature scale appropriately. To study the ef- 
fects of restricting the hopping to nearby sites only, we 
have studied a model where hopping is limited to sites 
within a radius r^- One would expect the parameter 
to decrease with density, so as to maintain a reasonable 
coordination number, z = (47r/3)nMn?'h- We use a phe- 
nomenological formula z = \2nMnl {nc + nMn), where ric 
is the Mn concentration for which the density of holes 
umi = pn-c corresponds to the metal-insulator transition 

1 /3 

nj^jjUB — 0.25. This formula has the proper asymptotic 
limits that z is proportional to n Mn at low density, satu- 
rates to z = 12 at large umu and at the Metal-Insulator 
Transition the average coordination number is z = 6, 
the same low coordination number as in the simple cu- 
bic lattice, and which is know to give a reasonable valug 
for the Metal-Insulator Transition for hydrogen centers.El 
We obtain Vh = 2.22aB for x = 0.0092 corresponding 
to an average z = 4.7, whereas for x — 0.05 we have 
Th = 1.59aB and z — 9.23. 

Another issue concerns the sign of the hopping inte- 
gral. In the model studied, the hopping integrals have 
been taken to have the same sign. However, the real sys- 
tem is heavily compensated due perhaps to As antisite 
defects, or to some other, as yet unknown, source. In any 
case, there will be potentials due to these compensating 
centers, which because of their opposite sign, may render 
some hopping elements of the other sign. In an effort 
to see the influence of such effects, we have also studied 
models where the hopping integrals have random sign. 

Figure |l^ shows the magnetization curves obtained for 
various models, as described in the figure caption. As 
can be expected, details of the hopping parameter lead 
to different Tc, consequently we show the data on a scaled 
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plot in terms of T/Tc- The crucial issue appears to be 
the density of states near the Fermi level (on the scale 
of the transition temperature Tc). Impurity bands that 
are broad have low Tc, while those that are relatively 
narrow, as might be expected for centers that have a 
strong central cell or on-site exchange energy, have high 
Tc. 

Despite the variation in Tc, we find a number of fea- 
tures that are generic to the models studied, as illustrated 
in Fig. |l| : 

(a) For each model t[r) studied, within the mean 
field approximation^ the Tc of the disordered system is 
higher than that of the corresponding ordered superlat- 
tice, while the magnetization at low temperatures ( below 
Tc/3 or so ) is lower for the disordered system. 

(b) The magnetization curves obtained with impurity 
bands from a set of positional disordered impurity cen- 
ters have an unusual shape (concave upwards), or at least 
significantly less concave downwards than the standard 
concave downward (convex upward) form of most uni- 
form magnets. 

(c) There is significant temperature dependence of 
M(T) at temperatures which are well below Tc, unlike 
in essentially uniform magnetic models where the mag- 
netization has reached its T = saturation value at Tc/2, 
or certainly by Tc/3. 

(d) Some of the unusual features of the Mn spin curves, 
which determine the bulk magnetization, are seen in the 
hole magnetization curves as well (see Fig.|l^ , but they 
are less pronounced. 

Other approximations made in our study involve ne- 
glecting a number of complications present in the experi- 
mental system such as: (i) the effect of the compensating 
centers, (ii) carrier-carrier interactions, and (iii) the va- 
lence band states. For (i), we have considered (subsection 
D) one of the effects, namely, if we include the random 
potential from these charged centers as an effective on- 
site random potential in our tight-binding Hamiltonian 
for the extreme cases (a) where the onsite potential is un- 
correlated, and (b) where it is modeled as being due to 
a close by As antisite defect. However, the longer range 
nature of the Coulomb potential in these systems with 
poor screening could also affect the hopping parameters, 
as we discussed in the preceding paragraphs. 

The states near the Fermi level for the range of con- 
centrations we consider are not strongly localized in the 
model studied, and even the actual system is fairly close 
to a metal-insulator transition (which we view as primar- 
ily driven by disorder, not electron correlation). Conse- 
quently, we believe that approximation (ii) is reasonable. 
An on-site interaction U (as in a Hubbard model), if 
treated in a mean-field approximation, would lead to in- 
creased Tc, since it aids in the splitting of the up and 
down spin bands. We show results for one such study 
that we carried out in Fig |l^ which confirms the above 
expectation. The change in Tc is about 35 % for the 
X = 0.01 case and 18% for the x = 0.05 case, for a value 
U = 1 Ry typical for hydrogenic centers at low densi- 



ties. However, this is likely an overestimate, since a fair 
fraction of the Mn impurity binding energy comes from 
short range potentials and exchange, and further, the ef- 
fective U near the metal-insulator transition is likely to 
be reduced by screening processes. 

Finally, we discuss the neglect of the valence band 
states. In the model we studied, we believe it is not im- 
portant, at least at low temperatures, where the shape 
of the magnetization curves is anomalous and where the 
mean field results should be at least qualitatively cor- 
rect. This is because the Fermi level lies ^ 300 meV 
above the valence band minimum, and so excitation to 
the valence band states would become dominant only at 
higher temperatures. It should be emphasized that the 
impurity states that we consider are derived from the 
host band (this would be the valence band in the case 
of Gai_a;Mn2;As). Consequently, if one wishes to include 
the rest of the band states of the host, they must be or- 
thogonalized to the impurity states; this has the effect of 
pushing the host band further from the impurity band, 
making the effect of host band states smaller than might 
be normally expected. In a more realistic model, how- 
ever, this remains an open question. 

IV. SUMMARY AND CONCLUSIONS 

In this study we analyzed a simplified model of III-V 
Diluted Magnetic Semiconductors, in which the charge 
carriers are restricted to a band formed from impurity 
orbitals, as in shallow doped semiconductors. We believe 
that this is a good starting point, because the high com- 
pensation present in these systems leads to carrier densi- 
ties which are not large enough to screen out the Coulomb 
interaction between the Mn ions and the charge carri- 
ers. This is clearly indicated by the presence of Metal- 
Insulator Transitions within the range of doping densities 
studied {x= 1-5% Mn). 

We started by analyzing an ordered superlattice case, 
with a homogeneous charge carrier distribution. This 
allows us to understand the rather unusual shape of the 
magnetization curves found especially for low values of p. 
We find good agreement with experimentally measured 
Tc with nominal parameters given in the literature. How- 
ever, we believe this is somewhat fortuitous, as it is well 
known in most magnetic models that mean-field analyses 
significantly overestimate the true Tc. Given the rather 
large number of free parameters, the experimental un- 
certainties about their exact values and the exponential 
dependencies on some of them in our model, it would 
always be possible to obtain good agreement with the 
experimentally measured Tc by changing the input pa- 
rameters within their error bars. 

We then show that positional disorder and on-site ran- 
dom interactions lead to significant changes in the shape 
of the magnetization and the critical temperature Tc, at 
least within a mean-field approximation. This enhance- 
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ment of the critical temperature by disorder has been 
confirmed, with differjeiit-techniques, in two recent stud- 
ies of the III-V DMS.EyO The critical temperatures ob- 
tained within our mean-field scheme for the same input 
parameters are larger than the ones measured experi- 
mentally by about 50% for the x = 0.05 sample, and 
about a factor of 3 for the x = 0.01 sample. While this 
may appear problematic, it is well known that mean- 
field schemes often significantly overestimate the true Tc, 
as we discuss later. It should be emphasized that the 
exchange coupling we have used, obtained from Bhat- 
tacharjee et alfl is lower than the values used by Mac- 
Donald et alJj or Millis and Das SarmaEHl for obtaining a 
similar Tc- The reason for this is that the impurity wave- 
functions that give rise to our impurity band are peaked 
at the Mn sites, and therefore provide a greater charge 
density at the Mn site than do host band wavefunctions. 

We have not performed detailed numerical fits to Tc 
because of a number of factors, such as (a) the compli- 
cated nature of the hole wavefunctions, (b) uncertainties 
about the nature of the compensation processes, (c) un- 
certainties about the precise description of the hopping 
integral, (d) the under-estimation of the effect of ther- 
mal fluctuations by the mean-field approximation, etc. 
One of the robust results that follows from our study 
is that the magnetization curves M (T) may vary con- 
siderably from the canonical concave downward (convex 
upward) form seen in practically all uniform magnetic 
models, independent of dimension or spin-components. 
Using different hopping parameters, we find that while 
Tc changes with the model chosen, the concave-upward 
form of M(T) over much of the temperature range below 
Tc is dependent mostly on the Mn spin concentration x 
and the carrier density px. The curves vary from essen- 
tially concave (convex upwards) functions (as reported in 
Ref. 1) for large relative charge carrier concentrations p, 
to almost linear dependence (as reported in Ref. 4), to 
very concave upward functions (as reported in Ref. 10) 
as p is decreased. Below T^ our calculated magnetization 
curves generically show a fast increase, followed by sat- 
uration and then fast increase again at much lower tem- 
perature, similar to the one reported in Ref. 19. Overall, 
we claim to find good qualitative agreement with the ex- 
perimental behavior concerning possible shapes of mag- 
netization curves, the metal-insulator transition etc. Our 
study suggests that by appropriate tuning of various pa- 
rameters, one may tailor the magnetic behavior M(H,T). 
in a manner not possible in simple uniform magnets. 
We believe that detailed information provided by experi- 
ments (using local probes such as ESR and NMR) would 
allow a clearer understanding of the nature of fcrromag- 
netism in these compounds. 

Returning to the issue of the magnitude of Tc, we be- 
lieve the most important correction to our mean-field re- 
sult is due to thermal (temporal) fluctuations, which are 
not included in a mean-fleld treatment. This can cause a 
substantial decrease in the critical temperature Tc- While 
typical renomalizations of Tc for uniform models range 



from tens of percent to factors of 2 or so, we expect them 
to be significantly larger for models with great spatial in- 
homogeneity, where percolation aspects may have consii. 
erable influence. Preliminary Monte Carlo simulationsEfl 
which include these fluctuations show that the the critical 
temperature is signiflcantly suppressed with respect to 
the mean-field value, especially for the lower concentra- 
tions X, where we found that the magnetization is some- 
what less robust. 

In addition to fluctuation effects, there are several in- 
teractions not included in our model which may lead to 
an overall decrease of Tc even within mean-field. One 
of these is the hole-hole Coulomb repulsion, which could 
prevent the accumulation of all the charge carriers in the 
high density regions, and would favor a more uniform 
spreading of the charge carriers over the entire sample. 
This would result in a smaller enhancement of Tc in the 
positionally disordered model vis-a-vis an ordered super- 
lattice of magnetic ions, and thus lead to a lowering of 
Tc- However, the magnitude of this effect may be small 
because the density of charge carriers is low. 

Another possible interaction not included in our sim- 
ple model, that may be more effective in preventing the 
concentration of charge carriers in the high-density re- 
gions, is direct Mn-Mn interactions. Mn-Mn interactions 
are expected to be AFM, as they are in II- VI DMS,c3 
while the charge carriers promote effective ferromagnetic 
Mn-Mn interactions. As a result, frustration is expected 
to appear, and to be most significant in the high-density 
regions where Mn-Mn interactions would be largest. In 
particular, nearest-neighbor Mn impurities may lock in a 
singlet state, and therefore not contribute at all to mag- 
netization. To first approximation, one may equate a 
system with such singlets to a system whose effective Mn 
concentration is smaller than the nominal one (the dif- 
ference being the concentration of singlets), and which is 
restricted to not having any nearest-neighbor Mn sites. 
In other words, it is as if these singlets become invisi- 
ble, as far as the magnetic properties of the system are 
concerned. Another phenomenon that may be responsi- 
ble for an effective homogenization of the sample is the 
creation of MnAs clusters, which is thought to be the 
reason behind the saturation of Tc for x > 0.05.a Al- 
though one can dope more Mn into the sample, an in- 
creased concentration x does not necessarily imply an 
increased number of magnetically and electrically active 
centers. In fact, some of the Mn impurities form small 
disordered six-fold coordinated centers with Asj— as ob- 
served in the closely related InMnAs compound£a Such 
centers are believed to give rise to n-type conductivity; 
in other words, they provide one possible compensation 
mechanism. These Mn impurities do not participate in 
the AFM exchange with the charge carriers, described 
above. One expects that the probability of appearance 
of such complexes is higher in the high-density Mn re- 
gions. Appearance of such complexes in the high-density 
regimes would lead to an effective decrease of the local 
active Mn density, and would effectively promote again 



12 



a more homogeneous ground-state, leading to a lower Tc. 
Finally, as the temperature increases towards Tc, one ex- 
pects that charge carriers may be excited to the valence 
band states (whose existence has been neglected in this 
model). This is likely to lead to a decrease of Tc, because 
the Bloch states would (a) lead to a more uniform dis- 
tribution of carriers, and (b) have less amplitude at the 
Mn sites than the impurity states. The magnitude (and 
importance) of this effect would, however, depend on the 
density of states of the impurity band and its separation 
from the host band. 

In conclusion, the nature of ferromagnetism in doped 
DMS is strongly affected by disorder. Surprisingly, 
we find that disorder actually results in a higher Tc. 
Whether this is due to the mean-field approximation, or 
is a more robust phenomenon, awaits results of Monte 
Carlo simulation studiesEj. Nevertheless, the versatility 
provided by a magnetization curve (and ensuing ther- 
modynamic properties) that has a tunable shape, makes 
DMS ferromagnetism a very interesting problem from a 
theoretical point of view. Adding to the richness are pos- 
sible effects of direct Mn-Mn interactions in concentrated 
systems (wMch lead to spin-glass behavior in undoped 
II- VI DMSE3), the existence of a ferromagnetic metal- 
insulator transition (unlike conventional doped semicon- 
ductors and amorphous alloys), and the likely unusual 
electron and spin transport characteristics because of dis- 
order. 
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FIG. 1. The average Mn spin Smii and the average spin 
per charge carrier = ^-Sh{i)/Nh = Suoie/p for dop- 
ing concentrations x = 0.01,0.02,0.03,0.04 and 0.05, and 
p = 10%. Due to their antiferromagnetic interaction, the 
two expectation values have opposite sign, with the Mn spin 
saturating at ^ and the charge carrier spin saturating at — ^ 
at low temperatures. The Mn spin curves have an uncharac- 
teristic shape, with inflection points as explained in the text. 
In the inset we plot the specific heat per Mn impurity as a 
function of temperature. The lower peak is due to Mn spin 
fluctuations, while the upper peak is due to charge carrier 
fluctuations. 
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FIG. 3. Comparison between experimentally measured 
critical temperatures Tc (full circles, from Ref. 10), and 
values obtained in the mean-field approximation for the or- 
dered Mn spin case. The experimental points fall between 
the lines corresponding to p = 5% and p — 10%, which is 
the estiipated range of charge carrier concentration for these 
crystals.Ej However, the Tc values are very strongly dependent 
on p (see Fig.^) as well as other parameters of the model. 
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FIG. 2. Evolution of the average Mn spin Smti and the 
average spin per charge carrier Sh for different charge carrier 
concentrations p — 5, 10, 25 and 40%. The Mn concentration 
is fixed at a; = 0.02. With large p, the shape of Smu becomes 
more like the usual convex upward function seen in typical 
ferromagnets. However, for low doping concentrations p the 
relatively few charge carriers can only fully magnetize all the 
Mn spins at very low temperatures. 
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FIG. 4. Left panel: Comparison between values for the 
average Mn spin Smu and the average charge carrier spin Sh 
obtained with biased initial conditions (full line) and random 
initial conditions (circles). For random initial conditions we 
actually plot IS'MnI > (full circles) and —\sh\ < (empty 
circles), since the two expectation values always have oppo- 
site sign, but the orientation is arbitrary. Results are for 
X = 0.00926, p = 10%. Right panel; the difference in energy 
per Mn spin between random and biased initial configura- 
tions, as a function of temperature. Full line is ksT /Nd- For 
details, see text. 
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FIG. 5. (a) Histogram (in arbitrary units) of the total den- 
sity of charge carriers at each Mn site logiQ(p(i)), for random 
Mn configuration with x = 0.00926 and p = 10%. The dot- 
ted line corresponds to ksT / J = 0.0006, while the full line 
corresponds to ksT/ J = 0.6. The vertical line indicates the 
position of the (5-function that describes the same histogram 
for an ordered Mn lattice. The two-peaked structure of the 
histogram for the disordered samples shows that in this case 
some Mn interact with holes much more strongly than the 
average, while some Mn spins interact with holes much more 
weakly than the average. Inset: Same, but for temperatures 
ksT/J = 0.6 (full line) and ksT/J = 1.5 (dotted line), (b) 
Comparison between the density of charge carrier histogram 
for X = 0.00926 (dotted line) and x = 0.05 (fuU line) systems. 
Both systems have p = 10% and fcsT/J = 0.0005. Vertical 
lines show the corresponding values for the ordered systems. 
The double-peaked structure is also visible in the higher den- 
sity sample. 
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FIG. 6. The average Mn spin Smu and average spin per 
hole Sh for doping concentration x = 0.00926 and p = 10%, 
In increasing order of Tc, the curves correspond to ordered, 
weakly disordered, moderately disordered and completely 
random distributions of Mn (see text). 
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FIG. 7. The average Mn spin Smu and the average spin per 
hole Sh for doping concentration x = 0.05 and p = 10% (thick 
lines) and 30% (thin lines) for typical random Mn distribu- 
tions (full lines) and simple cubic ordered Mn distributions 
(dashed lines). 
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FIG. 8. Magnetizations for 25 realizations of Mn posi- 
tional disorder, for x — 0.00926 (upper panel) and x = 0.05 
(lower panel). The long tails near Tc (especially for the lower 

density) arc duo to polarization of a few very dense clusters 
of Mn, and would be destroyed by thermal fluctuations. 
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FIG. 9. Average Inverse Participation Ratio as a function 
of energy for the a =1 (full symbols) and cr =t (empty sym- 
bols) sub-bands of completely random Mn distributions with 
p = 10%, X = 0.00926 and Na = 125 (circles), 512 (squares) 
and 1000 (triangles). The vertical line shows the Fermi en- 
ergy. Clearly, all the occupied states are localized. In the inset 
we plot the total Density of States for the two subbands, in 
arbitrary units. 
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FIG. 10. Average Inverse Participation Ratio as a function 
of energy for the a =1 subband of systems with strong dis- 
order (full symbols) and medium-disorder (empty symbols), 
for p = 10%, X = 0.00926 and Na = 125 (circles), 512 
(squares) and 1000 (triangles). The vertical line shows the 
Fermi energy of the system with strong disorder (full line) 
and medium disorder (dashed line). For these parameters, 
a medium-disordered sample is just above the MIT, while a 
strongly disordered one is just below the MIT. 



16 




Energy (meV) 

k^T/J 



FIG. 11. Average Inverse Participation Ratio as a func- 
tion of energy for the a =[ (full symbols) and a =] (empty 
symbols) subbands of completely random Mn distributions 
withp = 10%, X = 0.05 and = 548 (circles), 981 (squares) 
and 1600 (triangles). The vertical line shows the Fermi en- 
orgy, above the mobility edge. In the inset wo plot the total 
Density of States for the two subbands, in arbitrary units. 



FIG. 13. Comparison between magnetization curves in the 
absence of an on-site interaction u(i) = (dashed line), and 
with an on-site interaction due to As aiitisites, as explained in 
text (full line). The curves correspond to the same disordered 
positions for the Mn sites, and x = 0.00924 and p = 10%. 
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FIG. 12. Mn spin and charge carrier average magnetiza- 
tion as a function of temperature, for one random Mn config- 
uration corresponding to a; = 0.05 and p = 10%, and differ- 
ent values of the on-site disorder cut-off W/Ei, = 0, 0.5 and 
1. While on-site disorder does not affect Tc considerably, it 
does change the shape of the magnetization curves. 
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FIG. 14. Mn spin and charge carrier average magnetiza- 
tion as a function of temperature, for one random Mn configu- 
ration corresponding to x — 0.05 andp = 10%, for an external 
magnetic field H = 0,5 and lOT. The overall magnetization 
is significantly increased at lower temperatures, since the ex- 
ternal magnetic field polarizes all the Mn spins, not only the 
ones in the strongly interacting regions. 
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FIG. 15. Mil spill and charge carrier average magnetiza- 
tion as a function of T/Tc, for one random Mn configuration 
corresponding to a; = 0.00926 and p = 10%. The curves cor- 
respond to: (a) hopping allowed between all sites, same-sign 
hopping integral (Tc/J ~ 0-6); (b) hopping allowed only for 
sites with Vh = 2.22aB of each other, same-sign hopping in- 
tegral {Tc/J ^ 0.2); (c) hopping allowed between all sites, 
random-sign hopping integral [Tc/J ~ 0.12); (d) hopping 
allowed only for sites with rn ~ 2.22as of each other, ran- 
dom-sign hopping integral {Tc/J ~ 0.2). Sample averaging is 
required to obtain smooth curves. For more details, see text. 
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FIG. 16. Effect of an on-site interaction U — 1 Ry (treated 
within the mean-field approximation) on the magnetization 
curves corresponding to a; = 0.01 and p = 10% (dashed 
curve). For comparison, the curves obtained in the absence 
of an on-site interaction, for the same realization of disorder, 
are shown as well (full line). 
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